Modelling impurity-assisted chain creation in noble-metal break junctions 
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In this work we present the generalization of the model for chain formation in break-junctions, 
introduced by Thiess et al. [1], to zigzag transition-metal chains with s and p impurities. We apply 
this extended model to study the producibility trends for noble-metal chains with impurities, often 
present in break junction experiments, namely, Cu, Ag and Au chains with H, C, O and N adatoms. 
Providing the material-specific parameters for our model from systematic full-potential linearized 
augmented plane- wave first-principles calculations, we find that the presence of such impurities 
crucially affects the binding properties of the noble-metal chains. We reveal that both, the impurity- 
induced bond strengthening and the formation of zigzag bonds, can lead to a significantly enhanced 
probability for chain formation in break junctions. 
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I. INTRODUCTION 

Quantum transport in nanoscale and atomic-scale ob- 
jects has recently become one of the most researched 
areas in modern physics, owing to vast possible appli- 
cations in technology. On the most exciting side, elec- 
tronic properties of one-dimensional (ID) systems such as 
atomic chains and molecular wires often lead to fascinat- 
ing transport properties, which were predicted to occur 
theoretically and could be observed experimentally. Af- 
ter years of research, the material versatility of the class 
of ID systems still remains rather poor, however. This 
can be mainly attributed to the fact that to produce and 
sustain a free-standing, or even deposited chain proved 
to be rather difficult in experiment. One of the most 
promising experimental techniques developed over the 
last decade, which allows to produce ID structures and 
measure their transport properties, lies in investigating 
the mechanically controllable break-junctions (MCBJ). 2 

In a typical break junction experiment, two pieces of 
material being initially in contact are pulled apart while 
the conductance of the system is measured. During the 
process of pulling, several atoms may be extracted from 
the leads, forming a free-standing finite ID chain. 2 With 
this technique it has been possible to create monoatomic 
chains with the length of 5 to 10 atoms in Ir, 3 Pt 2 and 
Au 4 6 break junctions. The ability to create similarly 
long chains in break junctions of other metals has been 
greatly hindered so far. E.g. it has been demonstrated 
that the tendency for chain creation in lighter Ad and 
5d transition-metals (TMs) is significantly lower than in 
5d TMs. This trend was later explained based on the 
relativistic enhancement of mass in heavier elements with 



corresponding analogy to the reconstruction of the low- 
index surfaces of the 5d TMs. 

In the past years it became evident that it is possible 
to strengthen the bonds in a suspended chain and achieve 
higher degree of chain's producibility by adding external 
absorbates during the process of chain formation. This 
opens a new possibility for achieving desired versatility 
of produced in MCBJ experiments one-dimensional sys- 
tems with different properties. In a key experiment, 7 
Thijssen et al. observed such an enhanced tendency to- 
wards chain elongation for Au, Ag and Cu break junc- 
tions with oxygen present in the atmosphere. Moreover, 
the effect of other impurities, such as H, B, C, N and S, 
on the creation of Au chains has been considered, while 
stable and strong bonding between the Au and impu- 
rity atoms has been demonstrated. 8 10 Physically, such 
a situation can be explained based on a well-known ob- 
servation that low-coordinated noble metal (NM) atoms 
are chemically more reactive than in the bulk and thus 
their chains are expected to be even more reactive than 
nanopart icles . 1 1 ' 12 

Experimental advances in producing longer chains 
with the aid of impurities of various physical and chem- 
ical nature call for an appropriate theoretical descrip- 
tion. Taking into account the complexity of MCBJ ex- 
periments, the aim of the proper chain formation the- 
ory lies mainly in guiding future successful experiments 
and extending the material base for further studies. In 
this respect the predictive power of density functional 
theory (DFT) first-principles, or, ab initio, techniques is 
particularly valuable. Treating the complicated process 
of chain's stretching, elongation and eventual breaking 
upon pulling the electrodes apart in an atmosphere con- 
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taminated with various elements, is an extremely chal- 
lenging task even for the fastest first-principles methods. 
The main obstacle here is the necessity to describe si- 
multaneous restructuring of hundreds of atoms in real 
time. Due to computational limitations, ab initio stud- 
ies of impurity-assisted chain formation so far focused 
mainly on stability, mechanical and electronic properties 
of clean and contaminated infinite chains of transition 
metals. 13 16 This situation calls for a development of an 
effective model which is able to utilize the results of first- 
principles calculations of simplified systems for the pur- 
pose of describing the chain formation on a level, which 
allows to extend the material base for future MCBJ ex- 
periments with a reasonable computational effort. 

Recently, Thiess and co-authors introduced a model 
('CF model', or, 'chain formation model' in the follow- 
ing) able to provide a probability of chain's elongation 
in a break junction using only a limited number of pa- 
rameters which can be extracted from first-principles cal- 
culations of a perfect infinite chain and idealized leads. 
Using this model, the authors of ref. 1 analyzed the chem- 
ical trend of chain creation in break junctions of 3d, Ad 
and 5d transition metals, arriving at a qulitative agree- 
ment with experiments and unravelling the physics be- 
hind the observed trends. 1 An attractive feature of this 
approach lies in the ability to visualize the process of 
chain formation by referring to phase diagrams for prob- 
ability of elongation as a function of chain's stretching 
and the number of chain's atoms. Additionally, the ob- 
tained chemical trends can be understood based on a 
small number of element-specific parameters. Using this 
model, the tendency for chain elongation in linear Au-0 
chains has been predicted to be stronger than in pure Au 
chains, in accordance to experiments. 1 

In this work, we systematically apply the CF model 
for a detailed study of the trends in formation of Cu, Ag 
and Au chains in pure break junctions as well as Cu, Ag 
and Au break junctions with atmosphere contaminated 
by H, O, N, and C impurities. Additionally, we extend 
the CF model to the case of geometrically more complex 
zigzag arrangement of the atoms in the chain. In particu- 
lar, we find that some of the chains exhibit a finite proba- 
bility for elongation under significant tension in the linear 
regime as well as for the zigzag arrangement. While the 
probability of the chain's elongation in the zigzag regime 
is rather small, it has an advantage that it is more ro- 
bust with respect to breaking. Overall, we predict that 
incorporation of all considered impurities, even hydrogen, 
into the noble metal (NM) chains results in an increased 
probability for chain elongation when compared to the 
pure chains of NM elements. We also observe that al- 
though NM chains with s-type H impurities display wide 
regions of finite probability for elongation as a function of 
strain, this probability is much smaller than that for NM 
chains with p-type impurities. Remarkably, our calcula- 
tions predict a noticeably higher tendency towards chain 
elongation in Cu and Au contaminated chains, than in 
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Figure 1: (Color online) Schematic structure of (a) single- 
atom planar zigzag chains, (b) single-atom dimerized chains 
and (c) planar zigzag chains with impurities. The noble atoms 
are represented with large blue spheres, while smaller grey 
spheres stand for the impurity atoms. 

Ag chains. 

The paper is organized as follows. In the next sec- 
tion we provide the details of our DFT first-prinicples 
calculations. In Section III we analyze the energetics 
and structural properties of linear, zigzag and dimerized 
single- atom NM chains and NM chains with O, H, C, and 
N impurities. After this in Section IV we present the def- 
inition and the details of the model we use to estimate 
the tendency for chain elongation, as well as its gener- 
alization to the case of the zigzag geometry. Finally, in 
Section V we analyze the process of chain formation in 
break junctions of NMs, pure and with impurities, and 
compare the predicted trends. A summary and conclu- 
sions are given in Section VI. 



II. COMPUTATIONAL DETAILS 

For our density functional theory first-principles cal- 
culations we employed the full-potential linearized aug- 
mented plane- wave (FLAPW) method for ID systems, 17 
as implemented in the FLEUR code 18 within the gener- 
alized gradient approximation (GGA) to the exchange- 
correlation potential. The FLAPW basis functions were 
expanded up to /c max of 4.0 bohr -1 and we used 32 k 
points in one half of the ID Brillouin zone for our self- 
consistent calculations. The diameter of the cylindrical 
boundary between the interstitial and vacuum regions, 
Z} vac , and the fictitous in-plane lattice constant, D, were 
set to 8.3 bohr and 9.9 bohr, respectively. The muffin-tin 
radii have been set to 1.90 bohr for Au (with local orbitals 
for the 5s and bp states), 2.0 bohr for Ag and Cu, and 
1.0 bohr for the impurities. These values were chosen 
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Figure 2: (Color online) Upper panel: The minimal total 
energy e (per atom) as a function of d z for Au (left), Ag 
(middle) and Cu (right). The dashed red lines correspond to 
the total energy of the linear infinite single-atom chains, while 
the filled squares are the total energies of the dimerized chains. 
Lower panel: The evolution of the corresponding zigzag a- 
angles. 

to guarantee the required accuracy as well as to enable 
a comparison between the total energies of the zigzag 
chains with a wide range of the angle a, see Fig. 1. We 
have set the coordinate system such that the chain axis 
is aligned along the z axis, and considered a two- atoms 
unit cell to allow for zigzag and dimer arrangements. In 
Fig. 1 a sketch of the infinite planar zigzag and dimerized 
atomic chains is shown, where the definition of d z and d x 
distances, as well as the zigzag angle a are given. For 
our calculations we neglected the effect of the spin-orbit 
interaction, since its effect on the energetics of the chain 
formation was shown to be negligible. 19 ' 20 



III. ENERGETICS 

One of the purposes of our work is to understand the 
role of the geometric structure in chain formation taking 
into account energetics considerations. From this start- 
ing point we analyze the particular role of environmental 
impurities in the chain formation process. In order to 
achieve this goal, at first, we perform systematic total 
energy calculations as a function of the both interatomic 
parameters, d z and d x (c.f. Fig. 1). For each d z we per- 
formed a total energy calculation for at least five different 
d x values. On the basis of these calculations we inter- 
polated the two-dimensional energy profile in the d z -d x 
plane, from which we extract the effective minimum en- 
ergy curves as a function of d z , e(d z ) = min^ e(d z , d x ). 

Our results for the noble metal Cu, Ag and Au single- 
atom chains are summarized in Fig. 2, where the total 
energy and the zigzag a-angles are plotted as a function 



4 6 8 10 12 14 6 8 10 12 14 6 8 10 12 14 
d (a.u.) 

z v J 

Figure 3: (Color online) Upper panel: The minimal total 
energy (per unit cell) as a function of d z for Au (left), Ag 
(middle) and Cu (right). For better comparison, the energy 
curves were rigidly shifted with respect to each other. Lower 
panel: The evolution of the corresponding zigzag a-angles. 

of d z . The corresponding energies of the single-atom in- 
finite linear and dimerized chains are also presented for 
comparison. As can be seen from this figure, in all three 
cases the zigzag structures with large values of the a- 
angle are more stable than the linear configurations, in 
agreement with previous studies. 21 23 All energy profiles 
display a characteristic two- well structure, where the first 
minimum correspond to the a-angles in the vicinity of 
60°. In this configuration, each atom interacts strongly 
not only with its "chain" first nearest neighbors along 
the z-axis, but also with the second nearest neighbors. 
The second minimum in the total energy curves, which 
is clearly pronounced in particular for Au and Cu, oc- 
curs for the values of a at around 25°, corresponding 
to a zigzag structure with each atom interacting mainly 
with its two first nearest neighbors. While the zigzag 
angle falls off rapidly to zero for d z above 5.0 bohr for 
all chains, this decay is more abrupt for the case of Ag 
chains. As far as the dimerization is concerned, we find 
that the chains tend to dimerize when elongated above 
the d z of w 5.5 bohr, with a modest gain in energy as 
compared to the linear configuration for 5.5 < d z < 6.5 
bohr. As we shall see in the following, when the inter- 
atomic distance lies in this interval, the chains are already 
very prone to breaking, thus, we do not consider the ef- 
fect of dimerization on the chain formation considered in 
the following sections. 

Next, we follow the same procedure for the case with 
an impurity in the two-atom unit cell, as shown schemat- 
ically in Fig. 1(c). In Fig. 3 we present the evolution of 
the minimal total energy and corresponding zigzag angle 
with d z in Cu, Ag and Au chains with H, C, N, and O 
impurities (NM-X chains). Note, that the energy profiles 
are shifted in this figure with respect to each other along 
the y axis for better visibility. From Fig. 3 we observe 
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that the behaviour of the total energy is different as a 
function of d z for s- (H) and p-impurities (C, N, O). In 
the latter case the global minimum in energy occurs for 
the linear configuration (a = 0°) and the second local 
energy minimum can be seen also for the zigzag planar 
geometry with comparatively large angles a in the range 
of 30 — 40°. This is in contrast to the pure NM chains, for 
which two energy minima are also present, yet they both 
correspond to zigzag geometry with profound energetic 
preference to higher a-angles, see Fig. 2. On the other 
hand, for NM-H chains the zigzag structure exhibits the 
global and only energy minimum, although the values 
of the zigzag angle are rather small, namely, a « 12° 
(Au-H), a « 6° (Ag-H) and a w 10° (Cu-H). We can 
therefore conclude that the presence of impurities in NM 
break junctions leads to an effective straightening of the 
bonds in the chain with energetically favored configura- 
tions which are linear or exhibiting small zigzag angles. 

A common feature for all Ag-X chains, which can be 
seen in Fig. 3, is the smaller slope of the energy curves as 
a function of d z , when compared to corresponding Cu-X 
and Au-X cases. Such "softening" of the e(d z ) profile can 
be seen already for pure Ag chains, when compared to 
Au and Cu chains, especially in the vicinity of the global 
zigzag minimum, Fig. 2. The reason behind this phe- 
nomenon lies in the fact that the d states of Ag are situ- 
ated much lower in energy than the s-states, which solely 
determine the character of the electronic states in the 
vicinity of the Fermi energy (Ep). In case of Au and Cu, 
the d states are positioned higher in energy and overlap 
with the s states at the Ep, providing additional channel 
for hybridization and bonding between the atoms. Cor- 
respondingly, the bonding in the chains of Ag atoms is 
somewhat weaker than in those made of Au and Cu. Such 
weakening of bonding and softening of the energy profile 
for Ag in comparison to Au and Cu was also observed for 
linear infinite chains and in bulk, c.f. [24] and citations 
therein. 

In order to understand the difference in structural 
properties of NM chains with s- and p- impurities, we 
consider the density of states (DOS) of Au-0 [d z = 7.25 
bohr, a = 0°) and Au-H [d z = 6.0 bohr, a ~ 18°) chains 
in close to equilibrium cofigurations, presented in Fig. 4. 
In this plot, the DOS is additionally decomposed into 
atomic (Au or X) and orbital s, p and d contributions. 
By looking at the DOS of the Au-0 chain first, we observe 
that for the linear arrangement of the atoms in the chain 
the bonding of the atoms of Au and O is the strongest 
owing to the pronounced directionality of the oxygen's 
valence p-orbitals. In this case the overlap between the s 
and d z 2 orbitals of Au and p z orbital of O along the chain 
axis, as well as p x and p y oxygen orbitals with Au d states 
of the A3 symmetry (d zx and d yz orbitals) is the largest, 
which can be seen from the DOS (upper panel in Fig. 4). 
This leads to a pronounced minimum in the total energy 
for the linear regime. On the other hand, the hydrogen 
s orbital lacks directionality of the p orbitals of oxygen 
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Figure 4: (Color online) At omically- resolved density of states 
for Au-0 chain with d z = 7.25 bohr and a = 0° (top), and 
Au-H chain with d z — 6.0 bohr and a ~ 18° (bottom). The 
densities of states are additionally decomposed into the con- 
tributions coming from s, p and d states of different symmetry. 

atoms. In the linear case, brought into contact with Au 
atoms, the H s-state can hybridize only with the corre- 
sponding Au states of the same Ai symmetry, namely, 
6s and 5d z states, while the states of the A3 and A4 
(d xy and d x i_ y i) symmetry would remain non-bonding. 
Upon transition to a zigzag configuration, the cylindrical 
symmetry of the chain is broken and the hybridyzation 
between the hydrogen s and Au A3 orbitals occurs since 
it is not anymore forbidden by symmetry, as can be seen 
from the DOS in lower panel of Fig. 4. This additional 
channel for bonding in the chain energetically favors the 
zigzag arrangement of H and Au atoms, although the 
magnitude of the zigzag angle is limited by the small ra- 
dius of the Is hydrogen orbital. Noticeably, in both cases 
the A4 Au orbitals, which lie in the plane perpendicular 
to the chain axis, remain non-bonding, Fig. 4. 



IV. MODEL OF CHAIN FORMATION 

During the process of stretching and chain formation 
in a MCBJ experiment, the atoms are extracted from 
the leads upon pulling them apart, and form a chain 
consisting of several atoms. Owing to the complex re- 
laxations of hundreds of atoms during this process, its 
rigorous description from the first-principles is an ex- 
tremely challenging task. Moreover, the exact structure 
of the electrodes is normally unknown, which introduces 
a great deal of uncertainty in the initial configuration of 
the system. In the past, several ways to overcome these 
obstacles were taken, among which the most successful 
ones lie in employing the molecular dynamics technique 
to tackle the chain formation process. 13 ' 25 28 In this work 
we choose a different path, and extend a simple model 
of chain formation proposed by Thiess and co-workers, 1 
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which is based on general total energy arguments. The 
applicability of this model to description of the compli- 
cated MCBJ chain fomation process has been demon- 
strated for transition metals. 1 ' 19 

According to the CF model, we divide the system into 
two regions: the leads and the suspended chain. The 
electronic structure of the two parts is considered sep- 
arately from ab initio, thus neglecting their mutual in- 
fluence. This approximation becomes increasingly better 
when going to the limit of longer chains, which are in 
the focus of our study. The process of the chain forma- 
tion is monitored by following the change in the total 
energy of an atom that is transferred from the lead into 
the chain upon chain elongation. To estimate this to- 
tal energy change we assume that upon the transfer of 
the atom from the lead into the chain the only parameter 
which is changed in the system is the interatomic distance 
d between the atoms in the chain, while the length of the 
chain L and the structure of the leads remain the same. 
This is the so-called rapid reformation approximation 
(RRA). 1 ' 6 ' 29 ' 30 Moreover, we suppose that the transfer 
of a lead atom into the wire happens instantaneously at 
a certain time when the energy of the transfered atom 
and its coordinates change suddenly from their values in- 
side the lead to the corresponding values inside the chain. 
We also assume that between two consequent elongation 
events at times U and U+i the pulling of the electrodes 
apart results only in the smooth change of the chain's 
interatomic distance d. 

In Fig. 5 we show a sketch of the system under consid- 
eration, where the two figures for each of the situations 
correspond to the initial and final states of the process 
of the chain elongation. To analyze the producibility of 
the chain we compare the total energies of the initial and 
the final configurations in the process of pulling one lead 
atom into the chain, where L is the distance between the 
leads. A successful chain elongation will be translated 
into the following total energy relation, the criterion for 
producibility (P-criterion): 

(TV + l)e(d) > AE lead + (TV + 2)e(d) (1) 

where d = L/(N + 1) and d = L/(N + 2) are the chain 
interatomic distances before and after the elongation, re- 
spectively, and TV + 1 is the number of bonds in a chain 
of TV atoms. The binding energy of the suspended chain 
e(d) is defined from the binding energy of the infinite 
wire with interatomic distance d, Ew(d) = Ew(do) + e(d) 
(we always have two atoms in the unit cell), relative to 
the wire's cohesive energy at equilibrium interatomic dis- 
tance do j Ew(do), which implies e(d) > 0. 

The chain formation energy in eq. (1) is given by 
AEi ead = E w (d )-Ei ead , where E iead denotes the cohe- 
sive energy of an atom inside the lead. To calculate E\ ead 
we perform a bulk calculation of the total energy for each 
chemical element at its corresponding equilibrium lattice 
constant, Ebulk, and a bulk calculation of the total energy 



of the isolated atom, E^, to get AE bu i k = E^ - E bu i k . 
Taking into account the surface energies calculated by 
Skriver and Rosengaard, 31 AE sur f = E sur f — E ou i kl we 
defined E iead = AE bu i k - AE surf = E^ - E surf , where 
E sur f is calculated analogously to Ebuik for an atom at 
the surface of a certain crystallographic orientation. At 
finite temperatures, AEi ead corresponds to the difference 
between the chemical potentials of the lead and of the 
infinite wire. Typically, AEi ead > 0. If the P-criterion 
(eq. (1)) is satisfied the system tends to increase the num- 
ber of atoms in the chain by one. One can rewrite the 
P-criterion in the following way: 

£*(d,TV) = (TV + l)e(d) - AEiead + (TV + 2)e(d) (2) 

where E* (d, TV) is a function of the interatomic distance 
and the number of atoms in the chain. If J5*(d, TV) > 
the P-criterion is satisfied and the chain can increase by 
one atom for that particular pair of parameters (d, TV). 
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Figure 5: Sketch of the producibility (P) criterion in sus- 
pended single-atom chains (a) and in suspended chains with 
impurities (b). 

To apply the P-criterion to suspended chains it is use- 
ful to fit the chain binding energy e(d) with an analytical 
function. As the phase space of possible structural ar- 
rangements of the wires is two-dimensional, we calculated 
the binding energy for several d x and d z distances in the 
wire (see Figure 1), and fitted each energy curve with the 
analytical Morse potential, 1 interpolating the minimum 
energy curve from these data for all considered wires, as 
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explained in Section III. Thus, the binding energies par- 
ticipating in eqs. (1) and (2) are the energy curves shown 
in Figures 2 and 3. 

We define the string tension, or force F(d) as the slope 
of the binding energy with respect to the distance be- 
tween the atoms along the z-axis d = d z (see Fig. 5): 

5e(d z ) 



F(d) = F(d z ) 



(3) 



From eqs. (1) and (3) we can deduce that there are two 
main parameters which compete and which determine 
whether a chain can be increased by one atom. 1 These 
are (i) the energy cost of bringing an atom from the lead 
into the chain, AEi ea d, and (ii) the break force Fq, which 
is the maximum of the string tension occuring at the so- 
called inflection point d, that is, Fo = F(d). The inflec- 
tion point gives an estimate of how far an ideal infinite 
chain can be stretched until it breaks by the maximum 
break force due to long- wavelength perturbations. 1 Phys- 
ically, a successful chain elongation event will occur when 
the following happens: as the chain is stretched, the en- 
ergy of the system increases up to the point where a lead 
atom overcomes the chain formation barrier, AEi ea d, and 
enters the chain. This relaxes the distance in the chain 
from dtod and lowers the total energy of an atom in the 
wire, Ew(d) < Ew(d). The larger the slope of the total 
energy Ew(d), the more energy can be gained by relax- 
ing the chain from a distance d to d. Therefore, large 
values of Fq and small energy barriers AEi ea d will favor 
chain elongation. 

When the number of atoms in the chain N is very 
large, we can safely assume e(d) — e(d) w • Sd = 

F(d) • Sd, where 5d = d — d w jfe, and the P-criterion 
reads: 

F(d)> 



AE, 



lead 



•6(d) 



d 



/(d), 



(4) 



where f(d) is the generalized string tension, or the work 
done in drawing the chain with interatomic distance d 
out of the lead per unit chain length. This equation pro- 
vides the interval of distances in which the suspended 
chain can be producible. The analysis of the producibil- 
ity of the long wires and chains in terms of the general- 
ized string tension has been successfully performed in the 
past. 1 ' 19 ' 32 ' 33 It is easy to demonstrate that the lower and 
upper boundaries of the interval of distances where the 
chain is producible and eq. (4) is satisfied, correspond to 
the local minimum and maximum of f(d), respectively. 



V. IMPURITY-ASSISTED CHAIN FORMATION 
IN BREAK JUNCTIONS 



In this section we address the question of the impact of 
light s and p impurities H, N, C, and O on the producibil- 
ity of noble- metal chains in a MCBJ. A key towards un- 
derstanding the impurity assisted formation process is 



the additional degree of freedom of lateral relaxation in 
zigzag bonds, which can be incorporated into the crite- 
rion for producibility formulated in the previous section. 

Before analyzing the results of our extended model, it 
is important to discuss crucial energy scales determining 
the elongation process. While the extraction of atoms 
out of the leads can be modeled in line with previous 
studies, 1 ' 19 we have to consider additionally the supple- 
menting of the pure chain with the impurities. Here, 
we assume that the light impurities H, N, C, and O are 
present in the atmosphere and are not bound at the sur- 
face. This assumption is well justified considering that 
the concentration of such impurities is controlled by the 
partial pressure in experiments. Furthermore, a priori, 
the disociation energy of impurity molecules potentially 
acts as an additional barrier in the energetic balance 
given by relation (1). Here we refer to the results by 
Bahn et al. 12 who reported the chemisorption energy of 
Au chains on the order of —1.5 eV per O atom pair with 
respect to the energy of the free O2 molecule. This find- 
ing of highly preferential binding of impurity atoms to the 
chain allows for two important simplifications. Firstly, 
impurities have a strong tendency to cover all assisted 
bonds of the chain atoms. Hence, the resulting struc- 
ture is well described by a chain of alternating Au and 
O atoms. Secondly, we assume that the dissociation of 
impurity molecules and absorption to the metal bonds 
happens instantaneously after an additional metal atom 
has joined the chain. Thereby, we assume that the disso- 
ciation process happens independently of the chain for- 
mation and we do not account for it while considering 
the energetics of the chain formation. This leaves us with 
energy barriers introduced as AFi ea d in section IV, ac- 
counting for 0.86 eV, 0.75 eV and 1.16 eV for Au, Ag and 
Cu, respectively. 

Next, for an Au chain in the thermodynamic limit 
(N — >• 00) we introduce some important general con- 
cepts necessary for a detailed analysis of the results. As 
it is apparent from Fig. 6 and eq. (4) two distinct re- 
gions of producibility can arise: one in the domain of 
zigzag chains (angle a > 0) , and the second in the region 
of straight chains (angle a « 0). Accordingly, in both 
regimes the maximal sustainable force Fo is observed at 
two different values of the inflection point d. And while 
both parameters undergo large changes depending on the 
system, Fo and d of the straight chains have a clear phys- 
ical meaning and can provide important understanding 
of the most general trends of the chain formation. Hence 
— if not specified explicitly — Fq and d will refer to the 
straight chain domain from here on. 

The most obvious influence of considered impurities 
on the elongation process can be seen in the changes in 
the break force Fo when compared to the case of the pure 
chains. From Fig. 7, in which the break force for each of 
the studied cases is given, we first of all observe that 
in pure noble-metal chains a stronger binding and corre- 
spondingly larger values of Fq occur for Au and Cu, while 
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Figure 6: (Color online). Binding energy e(d), force F(cQ and 
generalized string tension f(d) for a pure zigzag Au chain. 
Both local maxima of F(d) are marked with open diamonds 
and will be refered to as break force Fo. Gray shaded areas 
are representing regions where f(d) > F(d) and the criterion 
for producibility is fulfilled. 



Ag chains have the weakest binding. This is a direct con- 
sequence of shallower energy profiles for Ag chains, dis- 
cussed in section III. The fact that Au chains exhibit the 
largest value of Fo among all NMs is in direct correspon- 
dence to experiments. 1 ' 2 When environmental impurities 
are added to the system, the bonds in all cases become 
stronger and the break forces larger owing to the pro- 
nounced covalent bonding in NM-X chains, as compared 
to metallic binding in pure chains. Noticeably, among 
all impurities H atoms cause smallest strengthening of 
the bonds, considering an absence of directionality of s 
orbitals, in contrast to the case of p impurities. On the 
other hand, adding p impurities to the NM chains leads to 
an approximately three- fold increase in the break force. 
Here, the most significant influence on Fo can be seen for 
NM chains with nitrogen impurities. 

Next, we perform a detailed analysis of the criterion of 
producibility for the NM and NM-X chains, as given by 
eq. (2). In Fig. 8 we plot the energy F* as a function of 
the interatomic distance in the chain, d, and the number 
of chain atoms, TV, both for pure NM chains, as well as 
NM chains with H, O, C and N impurities. The region of 
the (TV, d) phase space where F* > provides the param- 
eters of the chain configuration for which the elongation 
by one atom, or one NM-X pair of atoms, is possible. Ad- 
ditionally, the absolule value of the F* energy gives an 
estimate of how probable or how improbable the event 
of elongation is. Correspondingly, in Fig. 8 the absolute 
value of F* is given with appropriate color coding. The 
ability to estimate the F* value is extremely useful, since 
in experiments the energy barrier for an atom to enter 
the chain can fluctuate, depending on the geometry of 
the contacts. Therefore, high positive values of F* pro- 
vide higher credibility to the occurrence of a single event 



Figure 7: (Color online). Break force Fo for Au, Ag, and 
Cu pure and impurity- assisted chains. Along the x-axis "M" 
labels the pure chain, while "H", "C", "N", and "O" - chains 
with corresponding impurities. 

of chain elongation at a particular point in the phase 
diagram for which the P-criterion is satisfied. 

An important aspect of the chain creation process to 
consider is the stability of the suspended chains. The 
chain's stability in a break junction is a very compli- 
cated issue, since here the primary importance is due 
to the interplay of local tip geometry, phonon spectra 
of the finite chains, as well as temperature and applied 
bias. To model the chain's breaking process is therefore a 
highly non-trivial task, which is specific to each attempt 
of chain formation, and which cannot be undertaken in 
the present study. Here, we restrict ourselves to a simple 
but rather transparent way of quantifying the stability 
of the chains by referring to the value of their inflection 
point, d, i.e. the interatomic distance at which the max- 
imal force Fo is sustained. The value of which can be 
interpreted as the minimal interatomic distance at which 
the chain will be broken for a given TV, is presented for 
each of the cases by a thick vertical line in the phase 
diagram of Fig. 8. 

As we can see for pure NM chains in Fig. 8, the po- 
sition of the boundary of the P-region with respect to d 
is very similar for Cu, Ag and Au. On the other hand, 
for NM-X chains, the value of d varies in a compara- 
tively narrow interval of distance between 7.9 and 8.9 
bohr, Fig. 8. This means that the properties of consid- 
ered chains with respect to their stability are compara- 
tively constant. This observation is in accordance to a 
previous finding of ref. 1 in which the stability of linear 
suspended 5d transition metal chains has been analyzed 
within a better approximation. 

Let us first take a closer look at the phase diagrams 
for pure NM chains in Fig. 8. First we observe that all 
chains exhibit large areas of producibility, both for linear 
and zigzag geometries. This is a striking new finding, 
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Figure 8: Energy landscape E*(d, N) as a function of interatomic distance (d) and number of atoms in the chain (N) describing 
the chain producibility for pure and impurity-assisted noble metal chains, where E* is defined by the criterion for producibility, 
eq. (2). From left to right Au, Ag and Cu chains are displayed, from top to bottom, pure wires and H, C, N and O assisted 
chains. Thick white lines mark the E* = isolines indicating regions of producibility (E*(d,N) > 0). Additional thin isolines 
are drawn at E*(d, N) = 2 eV, 4 eV, 6 eV, and 8 eV in each panel. Note that regions which are outside of the interval [-8 eV, 
8 eV] are shown in the same color-coding as the minimal and maximal value of this interval. To rate the stability of chains, for 
each panel the inflection point d in the linear regime of the binding energy potential is depicted by thick blue vertical lines. 



which implies that in the general case of even more com- 
plex geometries of suspended wires more than one struc- 
turally stable phases favoring chain elongation can co- 
exists. This suggests that in experiment the co-existence 
of such phases can be used to enhance the elongation 



properties of the chains by switching from one phase to 
another via rapid pulling or pushing the electrodes apart. 
This technique can be particularly efficient, given that 
the stability of the phases and the respective elongation 
probabilities (proportional to E*) can depend strongly 
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on the details of the particular elongation attempt. 

Secondly, we observe that chains consisting of more 
than 10 atoms can be described very well in the ther- 
modynamic limit. Such length scale is rarely reached in 
experiment, 34 however, which underlines the importance 
of analyzing chain creation trends as a function of N. On 
this side, the description of the chain formation process 
at its initiation with one to three atoms in the chain is 
clearly improved by extending the CF model 1 to zigzag 
geometries: The strict assumption of linear bonds leads 
to an additional cost in energy for relaxations below the 
equilibrium interatomic distance. In our extended model, 
the chains relax at such short interatomic distances into 
a zigzag geometry with a comparable or even favorable 
energy contribution. Consequently, for Ag and Au the 
criterion for producibility can be fulfilled for smaller N. 
Generally speaking, this underlines the importance of 
structural versatility in a break junction experiment for 
chain elongation. Third, both linear and zigzag phases 
of producibility have largest values of E* energy for Au. 
The corresponding areas become less extended and ex- 
hibit somewhat smaller values of E* when going from Ag 
to Cu. This trend can be direclty related to the values 
Fo in Fig. 7 and its interplay with AE\ ea( ^. 

Going now to NM chains with H impurities (second 
panel in Fig. 8), we observe wide regions of producibility 
only in the linear regime for all metals, in accordance to 
small ground state zigzag angles in this case. And while 
the width of the P-region is roughly the same for Cu, Ag 
and Au, the probability of chain elongation, character- 
ized by the energy is noticeably larger for Au, than 
for Ag and Cu. In any case, the E* values, reaching as 
much as 5 eV for Au, are significantly larger than corre- 
sponding values for pure chains. Additionally, for all N 
the area of producibility can be entered upon increasing 
the interatomic distance in the chain without exceeding 
the inflection point. Overall, this leads to the prediction 
of a strongly enhanced tendency towards chain forma- 
tion and elongation in NM break juncction with H in the 
atmosphere. 

Finally, we turn to the case of NM chains with p im- 
purities C, O and N (three lower panels of Fig. 8). As 
we analyzed previously, adding p impurities to the NM 
chains results in a three-fold increase in the values of the 
break force (Fig. 7). This results in very wide regions 
of producibility for all chains with very high values of 
E* reaching over the threshold of 8 eV. These values are 
significantly larger than those for pure chains and even 
for chains with H adatoms. While the width of the P- 
region in the linear regime is similar for all NMs, the 
values of E* are somewhat smaller for Ag, than for Cu 
and Au, which can be attributed to the different break 
forces for these elements, see Fig. 7. Characteristically, 
for all cases the part of the P-region in the linear regime 
for which the interatomic distance lies below the inflec- 
tion point, is significant, thus ensuring stucturally stable 
elongation. 



As apparent from Fig. 3, contrary to the case of Cu 
and Au, Ag chains with C, N and O impurities exhibit 
two pronounced energy minimuma for the zigzag and lin- 
ear regimes, separated by a local energy maximum. In 
terms of the P-criterion this results in finite zigzag pro- 
ducibility regions for Ag, in addition to the P-region in 
the linear regime. And although the values of E* are 
modest for zigzag chains, the area in the phase space 
where the zigzag chains can be elongated is clearly vis- 
ible especially for C- and O-assisted Ag chains. As dis- 
cussed previously for pure chains, the second region of 
producibility can open the way to more efficient elon- 
gation. Despite the fact that shorter Ag chains can be 
grown in the linear regime only, for larger TV rapid push- 
ing of the tips towards each other might lead to a jump 
from the linear into the zigzag mode with further chain 
growth in the latter regime of higher stability. Clear indi- 
cation of the zigzag-like areas can be also seen for smaller 
d in AuC, AuO, CuN, and in particular CuC cases, al- 
though the energy E* is negative and reaches as much 
as —1 eV in these regions. Nevertheless, depending on 
the details of the chain growth and geometry of the tips, 
we believe that such regions can acquire positive energy 
E* and lead to longer chains in analogy to break junc- 
tions of Ag. Overall, we conclude that the tendency to 
elongation of NM chains with p impurities is enhanced 
when compared to pure NM chains and NM chains with 
hydrogen. 



VI. CONCLUSIONS 



We have presented the generalization of the model 
for formation of zigzag and impurity assisted chains in 
break junction experiments. On the basis of the extended 
model we investigated the growth of noble metal wires 
with H, O, C and N impurities, taking the input from 
first-principles calculations. Our most important obser- 
vation is that impurity assisted chain growth leads to 
a strongly enhanced tendency towards chain elongation, 
when compared to pure noble metal chains. Moreover, 
we reveal a distinct difference between the influence of 
the s and p type of impurities on chain creation, find- 
ing that p impurities lead to a higher chain elongation 
probability owing to a strong directional bonding. Im- 
portantly, we predict that the existence of more than 
one stable phases favoring chain elongation can occur in 
some cases, and suggest to employ this phenomenon for 
production of longer chains in experiment. By utilizing 
one of the most important advantages of our model we 
are able to translate our findings into phase diagrams for 
successful chain creation explicitly providing the values 
for the interatomic distances and the length of the chains 
for which the elongation can be achieved. 
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